Critical superfluid velocity in a trapped dipolar gas 
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We investigate the superfluid properties of a dipolar Bose-Einstein condensate (BEC) in a fully 
three-dimensional trap. Specifically, we estimate a superfluid critical velocity for this system by 
applying the Landau criterion to its discrete quasiparticle spectrum. We test this critical velocity 
by direct numerical simulation of condensate depletion as a blue-detuned laser moves through the 
condensate. In both cases, the presence of the roton in the spectrum serves to lower the critical 
velocity beyond a critical particle number. Since the shape of the dispersion, and hence the roton 
minimum, is tunable as a function of particle number, we thereby propose an experiment that can 
simultaneously measure the Landau critical velocity of a dipolar BEC and demonstrate the presence 
of the roton in this system. 



Liquid He was the first experimentally accessible sys- 
tem to exhibit dissipationless flow at low temperature, 
i.e., to demonstrate the existence of superfluidity in a 
quantum system. Landau famously explained this phe- 
nomenon by identifying a critical velocity vr, below which 
elementary excitations in the fluid could not be excited 
while conserving energy and momentum jjj. Because of 
this connection to the spectrum of elementary excita- 
tions, the Landau critical velocity can be expressed in 
terms of the fluid's dispersion relation oj(k) as 
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Remarkably, the Landau critical velocity vl does not 
coincide with the speed of sound in liquid helium, but 
is smaller due to the existence of an anomalously low- 
energy roton mode at wave vector k ~ A" 1 . This critical 
velocity was ultimately verified in experiments of ion drift 
velocity in liquid 4 He @- 

More recently, a new class of superfluids has been pro- 
duced in the form of Bose-Einstein condensates (BECs) 
of ultracold atomic gases. These gases have a distinct 
advantage over liquid helium in that they are dilute and 
hence easily characterized in terms of microscopic inter- 
actions. In particular, their critical velocity is nominally 
given by the speed of sound in the center of the gas, which 
can be easily calculated from the density and the s-wave 
scattering length of the constituent atoms. Early experi- 
ments at MIT sought to measure vl in a BEC of sodium 
atoms by stirring the condensate with a blue-detuned 
laser However, these experiments measured a crit- 

ical velocity for spinning off vortices rather than the true 
Landau critical velocity. This is a generic feature of such 
experiments in which the size of the object (in this case, 
the blue-detuned laser) is large compared to the healing 
length of the gas [5K3] . 

Still more recently, atomic BECs have been created 
whose constituent atoms possess magnetic dipole mo- 
ments large enough to influence the condensate [^, [l(| ■ 



These gases present a middle ground between atomic 
BECs and dense superfluid helium. Namely, the dipo- 
lar BEC (DBEC) is dilute enough to be understood in 
detail, yet its spectrum may exhibit roton features in pro- 
late traps, like those of liquid He 11]. The characteristic 
momentum of such a roton is set by the geometry of the 
trap in which it is held, whereas its energy is controlled 
by the density of dipoles, as well as the magnitude of the 
dipole moment [12] • Thus, by Eq. (p}, the Landau critical 
velocity is completely under the control of the experimen- 
talist. In contrast, vl in 4 He can be only weakly modi- 
fied by changing the pressure of the liquid [l3|. Thus, the 
DBEC provides an unprecedented opportunity to explore 
the fundamental relationship between the roton disper- 
sion and superfluidity. 

In this Letter we model an experiment on a DBEC 
similar to the MIT experiments. We consider a blue- 
detuned laser sweeping through a DBEC at a constant 
velocity, then compute the resulting condensate deple- 
tion due to the excitation of quasiparticles. We find an 
onset of depletion at a critical velocity that is near the 
Landau critical velocity at low densities. At higher densi- 
ties, where the roton determines vl, the critical velocity 
is a decreasing function of density, a behavior unique to a 
DBEC. Moreover, the simulations show a critical veloc- 
ity that is somewhat smaller than vj, at higher densities. 
We attribute this to the role that the roton plays in the 
mechanical stability of a DBEC. 

An ultracold, dilute DBEC containing N atoms is well- 
modeled within mean-field theory by the time-dependent 
non-local Gross-Pitaevskii equation (CPE), 
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where ^(r,^ is the condensate wave function, normal- 
ized to unity; r is the distance from the trap center; and 
U(r) = ^Muj 2 (p 2 + \ 2 z 2 ) is the cylindrically symmet- 
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ric harmonic trap potential with aspect ratio A = ui z /oj p 
where uj z and uj p are the axial and radial trap frequencies, 
respectively. The two-body interaction potential for po- 
larized dipoles with dipole moment d and zero scattering 
length is [14[ 



V(r - r') = d- 
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where 9 is the angle between r — r' and the polariza- 
tion axis. We choose the polarization axis to be the trap 
axis, z, so that the system is cylindrically symmetric. To 
characterize the strength of the dipole-dipole interaction 
(ddi) in a DBEC, we define the dimensionless quantity 

is the radial har- 



where a p = W h /Mo. 



monic oscillator length of the trap. The quantity D then 
characterizes either the density of the gas or the dipole 
moments of the atoms in the gas. 

We perturb this DBEC with a blue-detuned laser mov- 
ing at constant velocity v, which amounts to adding a 
potential 
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where a = 1 + (z/zq) 2 , zq = 7rwg/Ai as is the Raleigh 
length, wo is the beam waist of the laser, Ai as is the wave- 
length of the laser, y b (*) = ©(* — to)[v(t — to)] describes 
the motion of the laser in the y-direction and 6(t) is the 
Heaviside step function. This potential describes a laser 
that is stationary until t = to, at which time it moves to 
the edge with velocity v = vy. 

The effect of this blue-detuned laser on a DBEC is 
shown in Figure [T] for a DBEC with aspect ratio A = 
20, D = 124, and a laser with wo = 0Aa p and Uq = 
2Tiuj p where the chemical potential of the unperturbed 
condensate is fj, = 26.3huj p . We estimate the Landau 
critical velocity for this system to be vl ~ 1.5 a p uj p . For 
a laser velocity less than this (Fig.QJi), the condensate is 
completely unaffected whereas for a velocity larger than 
this (Fig. [Tp), quasiparticles are excited and the fluid 
would produce a net force on the moving laser. 

To determine the Landau critical velocity vl, we calcu- 
late the condensate's quasiparticle spectrum by solving 
the Bogoliubov de Gennes (BdG) equations [151 . Due to 
cylindrical symmetry of the system, the condensate plus 
BdG quasiparticles can be written as 

3 
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where u)j is the quasiparticle energy, m is the projec- 
tion of the quasiparticle momentum onto the z-axis and 
fi is the chemical potential of the ground state. Here, 




FIG. 1: The density profiles of a DBEC with D = 124 in a 
trap with aspect ratio A = 20 after a blue-detuned laser with 
axis z, beam waist Wq = 0.4a p , zo = 1.24a p and Uo = 2Tlu) p 
has traveled through the DBEC with velocity (a) v — 0.3 
a p uj p (b) v — 3.0 a p uj p . In (a), there are no visual excita- 
tions present in the system while in (b), excitations are clearly 
present, indicating the presence of a critical velocity for the 
system. The 1/e 2 contour of the laser is shown by the red 
dotted lines at the center of the condensates. 



ipa(p,z) is the stationary condensate wave function, i.e., 
the solution of Eq. ^ with time-dependence e _JMt , and 
is normalized to unity. The coefficients Cj (t) must be suf- 
ficiently small so that the BdG equations can be derived 
by linearizing the GPE about them. Their time depen- 
dence describes slowly varying quasiparticle occupations 
(compared to uj 1 ) in out-of-equilibrium states. 

In this formalism, the quasiparticles are characterized 
by their energies ojj and m quantum numbers. However, 
in order to apply the Landau criterion to this system, the 
quasiparticles must be characterized by a momentum, as 
well. To do this, we calculate the expectation value of 



the momentum, or (k p ) = ^J(k 2 ), of the quasiparticles. 

Using a Fourier-Hankel transform [HI , we transform the 
modes into momentum-space and compute the expecta- 
tion value of the linear momentum of the j th quasiparticle 
in momentum-space representation, 
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where we have time-averaged cross terms oc cos 2ujjt that 
oscillate on fast time scales [l6 ; ]. By associating these 
momenta to the excitation energies u>j, we determine a 
discrete dispersion relation for this system. The Landau 
criteria for superfluid critical velocity is derived by ap- 
plying conservation laws to translationally invariant flu- 
ids. Since the fully trapped system that we consider here 
is translationally variant, we apply the Landau criteria 
both to provide a hint as to where a critical velocity for 
quasiparticle excitations might be, and to test the appli- 
cation of this criterion to discrete systems. 

Figure [5] shows the discrete dispersion relations of a 
DBEC for various values of D. For D = (not shown), 
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FIG. 2: (Color online) The discrete BdG quasiparticle dis- 
persions for a DBEC in a trap with aspect ratio A = 20 for 
various values of D showing m = (black + sign), m = 1 
(teal squares) and m — 2 (pink circles) quasiparticles. As D 
is increased, the dispersion develops a phonon-like character 
at low momenta and a roton-like character at intermediate 
momenta. The slopes of the black dotted lines represent the 
corresponding Landau critical velocities for each D. 



teractions, and which corresponds to transverse sloshing 
of the condensate fl8| . The Kohn mode moves the con- 
densate's center of mass rather than exciting quasipar- 
ticles relative to the center of mass, which would imply 
the breaking of superfluidity in a translationally invari- 
ant system. We therefore ignore it here. In any event, we 
find that the occupation of the Kohn mode is very small 
compared to the total condensate depletion. 

We now compare vl as determined from the discrete 
dispersion relation with the onset of condensate deple- 
tion due to the laser having been moved through the 
DBEC. To quantify the breaking of superfluidity in the 
simulations, we calculate the depletion of the condensate 
by finding the quasiparticle occupations, or the number 
of particles that are excited out of the condensed state. 
In practice, this is achieved by calculating the ampli- 



tudes Cj (i) in Eq. ^ 
of the BdG modes 



via the orthogonality relations 
including their normalization 



Jdr k*(r>3'(r') - v*(r')v f {r')] = S jr , to give 
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where ^(r,^) is the numerical solution of the time- 
dependent GPE with the blue-detuned laser potential. 

The quasiparticle occupations are then given by rij (t) = 

n p uj p with (k p ) = y/n p + l/a p and \cj(t)\ 2 J dr'(|u 3 (r')| 2 + |v,-(r')| 2 ). In the simulations, the 

system evolves for a time T after the laser has completely 
left the system. We average the quasiparticle occupations 
for a time T after this, giving the average excited state 



the dispersion is given by the well known harmonic oscil- 
lator spectrum lo 



rip = 0, 1, 2, ... However, as D is increased, the spectrum 
changes to develop a phonon character at low-momenta 
and a roton character at intermediate momenta. In- 
deed, for D = 175.2, and more so for D = 230.0, there 
are some quasiparticles that branch off from the disper- 
sion towards lower energies and approach a momentum 
(k p ) ~ y/2Q/a p , corresponding to the characteristic roton 
wavelength A ro ton — 27ra z , where a z = yjTi/Mujz is the 



axial harmonic oscillator length [ll|, lUj ■ The modes with 
similar momenta but larger energy, on the upper branch 
of the dispersion, exist in lower-density regions of the con- 
densate while the quasiparticles on the roton branch exist 
in the high density center of the condensate. Note that 
Figure [5] includes only quasiparticles with m = 0, 1, 2. 

In each case, the Landau critical velocity vl is deter- 
mined according to Eq. (JXJ) as the slope of the shallow- 
est line through the origin that intersects a point on the 
dispersion curve; these lines are indicated in the figure. 
For smaller D, vl is determined by the low-momentum 
phonon-like modes where u> is linear in (k p ). By con- 
trast, for larger D, vl is determined by the low-lying 
roton mode and becomes a decreasing function of inter- 
action strength in contrast to a BEC with only contact 
interactions, where vl grows as the square root of scat- 
tering length. 

In evaluating vl from the discrete dispersion relation, 
we have ignored two excitations. One is the unphysical 
m = Goldstone mode. A second is the m = 1 Kohn 
mode, which has eigenvalue ui\ = Tiui p independent of in- 



occupations n,j — J Q dt'rij(t'). We find that T = 5 ui p 1 
is sufficient to converge these averages. 

Figure [3] illustrates the total quasiparticle occupation 
ntot = Y2j ™j as a function of laser velocity for vari- 
ous values of D using the laser parameters wq — 0.3a p , 
zo = 0.7a p and Uo = 0Ahui p . For each D, rh % stays very 
small until, at a certain critical velocity w cr it, it begins to 
increase significantly. Operationally, t> cr it is determined 
by the intersection of linear fits below and above u C rit- 
Well above u C riti the occupations decrease with velocity 
since the laser spends proportionally less time in the sys- 
tem as its velocity is increased. 

Notice that the overall depletion remains small with 
our weak laser. We have deliberately remained in the 
perturbative limit with our simulations to uncover the 
basic physics without the complications of large laser size. 
Additionally, we have checked that these lasers are not 
sufficient to excite vortex states in the DBEC. In prac- 
tice, larger condensate depletion would be obtained from 
a repeated back-and-forth stirring, as was done in the 
MIT experiments, or from a wider, stronger laser. While 
such a laser may spin off vortices in the condensate, thus 
defining a critical velocity smaller than Vl, the roton, for 
large enough D, would still determine the critical veloc- 
ity. 

Critical velocities determined from numerical Simula- 
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FIG. 3: (Color online) The occupations of the quasiparticles 
excited from a DBEC with aspect ratio A = 20 by a blue- 
detuned laser moving with velocity v (plotted on the horizon- 
tal axis) and with parameters too = 0.3a p , Uo = 0.4?kj p and 
20 = 0.7a p , for various values of D. At a critical v (indicated 
by the arrows), the occupations increase suddenly, indicating 
that the laser has excited quasiparticles in the system and 
superfluidity has been broken. 



FIG. 4: (Color online) The superfluid critical velocities v cr it 
for dissipation due to the excitation of quasiparticles in a 
DBEC as a function of D. The black dashed line represents 
the Landau critical velocity extracted from the discrete dis- 
persion relations of the system. The teal circles represent the 
results of numerical simulation for a laser with parameters 
Wo = 0.3a p , 20 = 0.7a p and Uo = OAfruip and the pink squares 
represent the results of numerical simulation for a laser with 
parameters wo = 0.4a p , 20 = 1.24a p and Uo = 2hw p . 



tions are presented in Figure|3]as a function of D. Results 
are shown for the comparatively weak (Uo = 0Ahu; p ) and 
strong (Uo — 2Huj p ) lasers. Also shown for comparison 
is vl (dashed line) as determined from the discrete dis- 
persion relations. At small D, the critical velocity grows 
slightly as the phonon modes stiffen and the speed of 
sound increases. This behavior is much like that of a 
BEC with purely contact interactions. 

At higher density, the critical velocity instead de- 
creases, due to the decreasing energy of the roton, and 
this is seen in both simulation and vl- The agreement 
is less perfect than in the phonon regime, however, with 
the simulated result coming in lower. This is because the 
roton, being the collapse mechanism for DBECs in traps 
with larger aspect ratios, softens with increasing conden- 
sate density. The presence of the laser in the DBEC 
serves to increase the density of the system, softening 
the roton and thus decreasing the critical velocity of the 
condensate, .just as a stationary laser leads a DBEC to 
instability [20j. For vanishingly small lasers, the criti- 
cal velocities extracted from numerical simulation show 
increasingly better agreement with vl- 

Finally, it is worthwhile to consider measurements of 
critical velocities in experimentally accessible DBECs, 
such as the 52 Cr system in Stuttgart @. Consider 52 Cr 
atoms whose scattering lengths have been tuned to zero 
in a trap with radial and axial frequencies uj p = lix x 100 
Hz and lu z — 27rx2000 Hz, respectively. This corresponds 
to a radial harmonic oscillator length of a p = 1.391 fim, 
particle numbers of N ~ 5701? and critical velocities 



in the range of 0.11 cm/s. These circumstances sug- 
gest that it may be plausible to observe the decline of 
the superfluid velocity with D for N > 8.5 x 10 4 52 Cr 
atoms, and hence to exhibit directly the roton 's influ- 
ence on superfluidity. This atom number corresponds 
to a maximum condensate density of n max ~ 9.5 x 10 14 
cm~ 3 , which, given the measured 3-body loss coefficient 
L3 = 2 x 10 -28 cm 6 /s [e|, should not produce significant 
losses over the time scales considered here. Additionally, 
we have checked that, for sufficiently large D, the roton 
serves to determine vl for 52 Cr DBECs with non-zero 
s-wave scattering lengths within the experimental uncer- 
tainty for 52 Cr, — 3ao < a s < 3ao [2l| . which is expected 
because these scattering lengths are sufficiently less than 
52 Cr's dipole length a dd ~ 15a [Hj]. 
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